Giant oscillations of energy levels in mesoscopic superconductors 
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The interplay of geometrical and Andreev quantization in mesoscopic superconductors leads to 
giant mesoscopic oscillations of energy levels as functions of the Fermi momentum and/or sample 
size. Quantization rules are formulated for closed quasiparticle trajectories in the presence of normal 
scattering at the sample boundaries. Two generic examples of mesoscopic systems are studied: (i) 
one dimensional Andreev states in a quantum box, (ii) a single vortex in a mesoscopic cylinder. 
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A normal cavity in a superconductor sample confines 
normal carriers due to their Andreev reflection from the 
walls formed by the superconductor order parameter. In 
the present Letter we show that Andreev levels in samples 
with sizes comparable to the coherence length exhibit 
giant mesoscopic oscillations as functions of the Fermi 
momentum kp and/or sample dimensions with an am- 
plitude that substantially exceeds the interlevel spacing 
they would have in bulk samples. For illustration, let us 
compare the effects of geometrical confinement for bound 
states in normal and superconducting systems. Consider 
one dimensional motion of a particle in a potential well 
of finite depth. The particle wave function oscillates as 
e ikx m c l ass i ca lly accessible- and decays exponentially in 
forbidden regions respectively. Placing the entire system 
into a quantum box of the size Lq larger than the width 
of the well, d, makes the wave function to vanish at the 
box boundaries. This is equivalent to a decrease in the 
effective width of the well and results in a slight modifica- 
tion of the bound states. Let us now take a sandwich- like 
structure of a total thickness Lq where the normal (N) 
slab of a thickness d is confined between two supercon- 
ducting (S) layers with a certain order parameter phase 
difference 4> between them (see Fig.^l, and consider the 
effect of the geometrical confinement on Andreev states 
with energies e below the superconducting gap A. As 
before, the particle/hole wave functions in the N region 
oscillate with k — kp ± e/hvp, where vp is the Fermi ve- 
locity. However, in contract to the previous example, the 
wave functions oscillate further into the S layers with the 
period of 2n /kp and with the amplitude slowly decay- 
ing on the scale of the superconductor coherence length 
£o = tiVp/A. Only when the particle in the box is exactly 
in the geometrical resonance, kpLo — im, i.e., when one 
of its normal-state energy levels coincides with the Fermi 
level, the Andreev states do not feel the external bound- 
aries, see Fig. ^ When the particle is out of resonance, 
sm(kpLo) ~ 1, to satisfy the zero boundary conditions 
an adjustment of either the wave vector in the N region 



FIG. 1: Andreev states in an SNS structure placed in a quan- 
tum box of length Lq. Solid (1) and dashed (2) lines show 
the positions of the box boundaries for resonance and off- 
resonance situations, respectively. 



by 6k ~ 5e/Hvp ~ 1/d (for d > £o) or of the phase shift 
between the electrons and holes by <5arccos(e/A) ~ 1 
(for d < £o) is needed. Thus, the deviations in energy 
are large and may compare to A for d ~ £o . 

We thus see that the interplay between geometrical 
and Andreev quantization results in giant oscillations of 
the Andreev levels as functions of kpLo with an ampli- 
tude of the order of A. The amplitude decreases expo- 
nentially as a function of the distance from the external 
boundary to the Andreev turning point of a particle at 
the NS interface. This oscillatory phenomenon can be 
viewed as a generalization of geometrical effects caused 
by the presence of impurity atoms in vortex cores Q 
and Tomasch oscillations in films and finite-size type- 
II superconductors 0, Q but with the enormously am- 
plified magnitude. The mesoscopic oscillations strongly 
affect both thermodynamic and transport properties of 
nanoscale superconductors that are the focus of current 
experimental and theoretical research (see 0- [J Q an d 
references therein). Interference effects of similar origin 
have been previously studied for hybrid NS systems with 
multiple semi-transparent insulator potential barriers [(| 
(see also Q for a review). 

Closed trajectories.- The one dimensional physics in 
a quantum box is related to the concept of closed tra- 
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jectories whose impact on the validity of the quasiclas- 
sical description of superconductivity is discussed in 0- 
Making use of this concept enables one to generalize the 
above picture to higher dimensions. A standard semiclas- 
sical approach to superconductors is formulated for the 
quasiparticle motion along the beams VS" where S is the 
normal-state eikonal, |V5| = kp. The wave function of 
fcrmionic excitations has two components in the particle- 
hole space W = (u,v) and can be written as ^ = -0e lS , 
where if) is the envelope function varying slowly over the 
Fermi wavelength. Generally, the phase S gained along 
open trajectories does not affect excitation characteristics 
such as energy levels or density of states (DOS). However, 
the closed trajectories formed due to boundaries, impu- 
rity potentials, and/or in an applied magnetic field, do 
influence the energy spectrum. 

Consider a closed trajectory of the length L. The wave 
function ^ should be single valued, which gives 

i>(s) = $(s + L)e lS ^ , (1) 

where S(L) = § VS 1 • dr = kpL is independent of the 
arc length s along the trajectory. Equation Q suggests 
that the initial problem is equivalent to the problem of 
an unbounded motion of a particle in the periodic gap 
potential A(s) = A(s + L) for the proper choice of the 
quasimomentum q (see below). In the latter case the 
wave function satisfies the Bloch theorem 

Ms + L) = e^ q (s) , (2) 

while the energy is a periodic function of q: e(q + 
2ir/L) = e(q). Comparing Eqs. and (J2J) we find 
q = [2ttM — S(L)]/L, where M is a large integer cho- 
sen such that q belongs to the first Brillouin zone. Con- 
sider two examples: (i) States with e > A = const will 
have the spectrum e 2 = A 2 + H 2 v f (2ttM' / L — kp) 2 where 
M' = M — N, N being the number of the energy band. 
This spectrum results in Tomasch oscillations [2, HI • (h) 
For sub-gap states the correspondence between the spec- 
trum in a mesoscopic superconductor and in a bulk sam- 
ple can be easily established for L 3> £o- Let e' ' be 
an Andreev bound state for L — > oo. The tight binding 
approximation in the equivalent periodic problem yields 

cos a + C) , (3) 

where the band width 5 ~ Ae~ AL is proportional to the 
exponential overlap of the decaying functions for the sub- 
gap states, A(e) ~ C^ 1 , and a = qL + (3. The parame- 
ters (3 and C depend on the particular problem. Since 
cos a — cos(kpL—f]), the energy level oscillates rapidly as 
a function of kpL. The amplitude of oscillations can well 
exceed the value of e^ ^ itself, provided the loop length L 
is not much larger than £o- 

The relative contribution of such oscillations to bulk 
properties depends on the relative weight of closed tra- 
jectories allowed by the particular sample geometry. Be- 
low we focus on two problems where these contributions 



are critical: (i) one dimensional (ID) Andreev states in a 
quantum box (this problem has been discussed briefly in 
the introduction) , and (ii) energy states in a vortex core 
placed in a clean mesoscopic cylinder of a finite radius. 
In both cases the dimensions of the system are assumed 
comparable to £o- 

ID Andreev bound states in a quantum box. - Con- 
sider a quantum point contact which is transparent to 
a few modes N c passing from one superconducting lead 
to another. We further assume that these modes are lo- 
calized within the device by specular reflections at the 
boundaries of the leads which are separated by a dis- 
tance Lq. The leads contain also large numbers of other 
modes Ni ea( i 3> N c , which take part in the superconduct- 
ing pairing. The leads can be connected, via some of the 
modes Ni eac i, to an external superconducting circuit to 
control the phase difference 4> between them. A possible 
realization of this device is an adiabatic constriction of 
the type discussed in ||. For the modes N c that pass 
through the constriction but are confined within the box 
of the size Lq, the quantum mechanical problem corre- 
sponds exactly to that in Fig.Q]for the limit d -C £o- 

We assume a step-like gap potential A(x) = 
A e ! sl s n ( a: ) 4>/ 2 . The confinement couples the states with 
opposite momenta and creates a closed trajectory loop 
of the length L = 2Lq. In the S region one has the 
waves e l ^ ±x and e~ l ^ ±x , where q± — k x ± iX, A = 
y Ag — e 2 /hv x , k x and v x are the particle momentum 
and velocity projections on the x axis. We choose here 
|e| < Ao, however, the same expressions hold also for 
|e| > Ao with an imaginary A. Matching these wave 
functions yields the dispersion relation 

e 2 = A 2 [l-Tsin 2 (0/2)] . (4) 

The transmission coefficient T = (1 + A)^ 1 , and A — 
sin 2 (fc :E Lo)/sinh 2 (ALo). For L > £o and e < A , Eq. (@} 
looks like Eq. in accordance with the general argu- 
ments above. Equation J3J has a familiar form Q of the 
spectrum for a contact with the double barrier of strength 
A. It describes mesoscopic fluctuations and accounts for 
the resonance transmission at sm(k x Lq) = 0. However, 
the effect of spectrum modifications is not restricted to 
mere renormalization of the barrier strength. Qualita- 
tively new features appear due to the energy dependence 
of the transmission coefficient T. Because of the geo- 
metrical quantization, the spectrum for |e| > Ao is no 
longer a continuum and cannot be separated from the 
phase-dependent sub-gap state; instead, we obtain a set 
of (^-dependent discrete levels in the entire energy range. 
For a short box Lq < £q, the lowest energy level is 

e 2 ((f>) = A 2 , [cos 2 (0/2) + (^ x /A i ) 2 sin 2 (fc ;CJ Lo)] ■ 

For |e| 3> Ao, the levels transform into the ^-independent 
spectrum in a normal-metal box: (hk x ±e n / v x )Lq = nnh, 
where n is an integer. Each ^-dependent level provides 
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oscillatory contribution to the supercurrent 0: 

/n(0) = _^%M tanh !!iM ; (5) 



where only e„ > are taken. Note that the supercurrent 
in Eq. JJJ is transported by the modes which, in the 
normal state, are localized and do not carry current. In 
the superconducting state, however, the current along 
the localized modes N c appears due to the conversion 
over distances ~ min{Lo,£o} of the supercurrent from 
delocalized modes Ni ea d- The analysis of Eq. (@J shows 
that the strongest <j> dependence is realized for the lowest 
energy state which thus dominates the current. 

Vortex core states in a mesoscopic superconductor. - 
We consider now low energy core states in a single vortex 
introduced into the center of a mesoscopic cylinder of 
a radius R > £o with the quasiparticle mean free path 
< > i?. Bogoliubov-de Gennes (BdG) equations read 
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where A = \A(r)\e t,Jz ^iay, <r, are Pauli matrices, and 
r, (j), z are cylindrical coordinates with the z axis parallel 
to the cylinder axis. The gap | A(r) saturates at Ao far 
from the vortex axis. The vector potential here is neg- 
ligible for an extreme type-II superconductor. We look 
for a solution *P = e l<Tz ^/ 2+v ^W with a given half-integer 
angular momentum \x. If the superconducting cylinder 
is surrounded by an insulator, the boundary condition 
requires U(R, z) = 0. 

We find the energy spectrum both analytically and nu- 
merically. For numerical computations we use a ma- 
trix representation of the BdG operator in the basis 
of the eigenstates with a given momentum k z along 
the z axis for a normal-metal cylinder of the radius 
R. We truncate the infinite matrix keeping the num- 
ber of eigenstates larger than the number of propagat- 
ing modes in the normal-metal waveguide. The obtained 
fc z -dependent matrix is diagonalized yielding the energy 
spectrum for a vortex. We approximate the gap by 
|A(r)| = A r/-\/r 2 + £ 2 , choosing = £ without the 
loss of generality. The calculated spectra do not depend 
qualitatively on the exact shape of |A(r)|. Shown in 
Fig. [21 are typical energy spectra calculated for a real- 
istic material parameter Ao/ Ep — 0.01. 

The analytical description is based on a standard qua- 
siclassical scheme [l£j modified to take account of the 
proper phases [lll | of radial waves for small values of fi: 



U = e lk * z H (1) 



(k r r) w^' 
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(k r r) w 



(-) 



where H, 



(1,2) 



arc the Hankcl functions, k; 



+ k 2 = k 2 



F ' 



and w = (w\, W2) are slow functions of r. Solving equa- 
tions for the envelopes w^' one can construct functions 




FIG. 2: The energy spectra (solid lines) for a vortex in a 
mesoscopic cylinder with (a) R/£o = 3.5 and (b) R/£o = 4.0 
as functions of k z . The CdGM energy spectra are shown by 
the corresponding dashed lines. 



w > and w < decaying at different ends of the trajec- 
tory passing by the vortex [12J. The wave function in 
a cylinder is the superposition w = A > w > + A < w < . 
The boundary condition at r = R couples the incom- 
ing and outgoing waves which leads to formation of a 
closed trajectory loop for one-dimensional motion along 
r. The solvability condition of the two linear homoge- 
neous equations for two constants A < and A < gives the 
bound state energy e(q) = e Al (/c r ) in the form of Eq. @ 
where S = A /(A cosh[2K(R)]), a = 2k r R + tt/2, 

f^ -* — = —uj(k r )fx is the Caroli-de 

for a vortex in a 



(7 = 0, and = e w = 
Gennes-Matricon (CdGM) energy 
bulk superconductor, 



uj(k r ) 



2mA 
JfkJX 



(|A(r)|/r)e~ 2K(r) dr , 



K(r) = 



\A(r')\dr', A 
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These analytical expressions are in a very good agree- 
ment with our numerical results. The last term in Eq. @ 
describes the mesoscopic level fluctuations. Their ampli- 
tude 8(k r ) is much larger than the CdGM level spacing 
u>{k r ) ~ Al/Ep if R is not exceedingly larger than £ . 
Different levels can cross each other because they belong 
to different angular momenta /1. The fluctuating levels 
also can cross zero for not very high fi. The amplitude 
5(k r ) decreases as k z approaches kp since cosh[2K(R)] 
in it grows exponentially as k r decreases. Note that in 
a layered two-dimensional superconductor the amplitude 
S(kp) is constant corresponding to k z — 0. For later use, 
we define the critical radius R c for which the maximum 
amplitude of oscillations S(kp) is equal to half of the max- 
imum distance between the CdGM states u>o = u>(kp). 
According to Fig. Efb) the radius R — 4£o is slightly 
larger than R c for A Q /Ep = 0.01. 

Density of vortex core states - Consider a sample ra- 
dius R < R c such that S(kp) ujq. According to 
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FIG. 3: Solid lines: The DOS for (a) R/£ = 2.5 and (b) 
R/Zo = 4. The CdGM DOS is shown by dashed line. 



Fig.^a), we divide the range < k r < kp into the region 
of large oscillations k* < k r < kp, where 5(k r ) > u>(k r ), 
and the region < k r < k* , where levels are smooth func- 
tions of k r close to the CdGM form (k r ). The total 
DOS is a sum v = v\ + 1*2 of contributions from the both 
regions. In the region k* < k r < kp, a large number of 
levels with different p cross the constant energy line many 
times so that the continuous approximation is appropri- 
ate. Therefore, V\ is equal to the CdGM DOS averaged 
over a large number of angular momentum cigcnstates: 



vi(K) = - 



k* Uj(k r )y/k 2 F - fc^ 



(6) 



(per spin projection). It provides a background zero- 
energy DOS due to mesoscopic level fluctuations. The 
term //2(e) is a sum of peaks positioned at e = e^\k*) 
with the energy period tj(k*) larger than the CdGM mini- 
gap wo- In the limit e ^> uj(k*) the term V2 (e) saturates 
at the value determined by the same integral as in Eq. © 
but taken within the limits to k*. Therefore, the total 
DOS saturates at the averaged CdGM DOS u = ^i(O). 

If R > R c but the amplitude of oscillations is still 
comparable to ujq, there are many roots of £p(fc r ) = e 
for a given p, such that the effective minigap is e m ; n = 
ojq/2 — S(kp). It vanishes for R = R c . 

To simplify our equations we approximate the CdGM 
interlevel spacing as uj(k r ) » (kp/k r )u>o while K(R) w 
(R/£o)(k F /k r ) with Awl. We find k*/kp = R/R c and 
Rc = (Co/2) In (Ao/wo). The period of DOS oscillations 
is thus u>(k*) = luqRc/R- The background DOS is 



V\ = Vq 



1 - 



(2/vr) fare 



csmp — py 1 — p* 



where Vq — kp/AujQ and p = Rj R c . It vanishes for p = 1. 

We calculated the DOS numerically using the obtained 
analytical expressions for the energy spectrum. To ex- 
clude a large number of van Hove singularities the DOS 
was averaged over a small energy interval Se = O.Iluq. 
The results for R < R c and R « R c shown in Fig. |3fa) 



and (b), respectively, are in good agreement with the 
above analytical estimates. In particular, the period of 
oscillations in Fig.^a) is approximately 1.5 times larger 
than the period u>o for the CdGM DOS in the bulk, which 
agrees with the value of the cylinder radius R « R c /1.5. 
At the same time, the period in Fig.^h) almost coincides 
with luq; in addition, the minigap here vanishes. These 
features well correspond to R being close to R c . 

To summarize, we predict a profound effect of geomet- 
rical quantization on Andreev states in mesoscopic super- 
conductors, which exhibit giant oscillations as functions 
of the particle momentum and the sample size. In partic- 
ular, the geometrical quantization results in appearance 
of zero energy modes for vortex core states. We discussed 
the case of ideal sample surfaces, however one expects all 
the essential conclusions to hold for atomically smooth 
sample surfaces as well. The spectrum oscillations can 
be observed by scanning tunnelling spectroscopy with 
high energy resolution and by transport measurements 
in weak links. 
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